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We analyze the stationary spin-dependent transport through a single-molecule magnet weakly 
coupled to external ferromagnetic leads. Using the real-time diagrammatic technique, we calculate 
the sequential and cotunneling contributions to current, tunnel magnetoresistance and Fano factor 
in both linear and nonlinear response regimes. We show that the effects of cotunneling are predomi- 
nantly visible in the blockade regime and lead to enhancement of tunnel magnetoresistance (TMR) 
above the Julliere value, which is accompanied with super-Poissonian shot noise due to bunching of 
inelastic cotunneling processes through different virtual spin states of the molecule. The effects of 
external magnetic field and the role of type and strength of exchange interaction between the LUMO 
level and the molecule's spin are also considered. When the exchange coupling is ferromagnetic, we 
find an enhanced TMR, while in the case of antiferromagnetic coupling we predict a large negative 
TMR effect. 

PACS numbers: 72.25.-b, 75.50.Xx, 85.75.-d 



I. INTRODUCTION 

Owing to recent advances in experimental techniques, 
it is now possible to study transport properties of individ- 
ual nanoscale objects, like quantum dotsji nanotubes^— 
and other molecules^— Investigation of electron trans- 
port through molecules is stimulated by the prospect of 
a new generation of molecule-based electronic and spin- 
tronic devices. It turns out that owing to their unique 
optical, magnetic and mechanical properties, molecules 
are ideal candidates for constructing novel hybrid de- 
vices of functionality which would be rather hardly ac- 
cessible in the case of conventional silicon-based elec- 



tronic systems 
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For instance, one interesting feature 



of nanomolecular systems, which does not have counter- 
part in bulk materials, concerns the interplay between the 
quantized electronic and mechanical degrees of freedom,- 
In this paper we deal with one specific class of 
molecules which possess an intrinsic magnetic moment, 
referred to as single- molecule magnets (SMMs) J^r— 
Such molecules are characterized by a significant Ising- 
like magnetic anisotropy and a high spin number S, which 
give rise to an energy barrier that the molecule has to 
overcome to reverse its spin orientation. At higher tem- 
peratures, the SMM's spin can freely rotate, whereas be- 
low a certain temperature it becomes trapped in one of 
two metastable orientations. Since magnetic bistability is 
one of the key properties to be utilized in information pro- 
cessing technologies, SMMs have attracted much atten- 
tion and a great deal of effort was undertaken to measure 
electronic transport through a SMMiiS - — The experi- 
ments carried out to date have concerned only the case 
of SMMs coupled to nonmagnetic electrodes. However, it 
has been suggested recently that spin-polarized currents 
(when the leads are ferromagnetic, for instance) can be 



used to manipulate the magnetic state of a SMM<22r— 
Such a current-induced magnetic switching (CIMS) of a 
SMM takes place as a consequence of the angular mo- 
mentum transfer between the molecule and conduction 
electrons. 

When considering coupling strength between the 
molecule and external leads, one can generally distin- 
guish between weak and strong coupling regimes. In the 
latter case, i.e. when resistance of the contact between 
the molecule and electrodes becomes smaller than the 
quantum resistance, the electronic correlations may lead 
to formation of the Kondo effector— These correlations 
result in a screening of the SMM's spin by conduction 
electrons of the leads, giving rise to a peak in the den- 
sity of states and full transparency through the molecule. 
On the other hand, in the weak coupling regime, the 
Coulomb correlations lead to blockade phenomenal For 
voltages lower than a certain threshold value, sequen- 
tial tunneling processes through the molecule are then 
exponentially suppressed due to Coulomb correlations 
and/or size quantization. However, once the bias volt- 
age exceeds the threshold value, the electrons can tunnel 
one-by-one through the molecule. The latter regime is 
known as the sequential tunneling regime, and the for- 
mer one is often referred to as the Coulomb blockade 
or cotunneling regime i 35 ' 36 It should be noted, however, 
that although the sequential processes are suppressed in 
the Coulomb blockade regime, current still can flow due 
to second- and higher-order tunneling processes, which 
involve correlated tunneling through virtual states of 
the molecule. Furthermore, although higher-order pro- 
cesses play a substantial role mainly in the cotunnel- 
ing regime, they remain active in the whole range of 
transport voltages, especially on resonance, leading to 
renormalization of the molecule levels and smearing of 
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the Coulomb steps^l Therefore a suitable theoretical 
method should be used to properly investigate trans- 
port through molecules in the regime where both the 
sequential and cotunneling processes coexist and deter- 
mine transport properties. The existing analytical studies 
of electronic transport through SMMs in the weak cou- 
pling regime were based on the standard perturbation 
approach^ ~ 28 i 38 i 39 and they dealt separately either with 
the sequential or cotunneling regime, with one attempt 
of combining them4£ Nevertheless, to properly take into 
account the nonequilibrium many-body effects such as 
for example on-resonance level renormalization or level 
splitting due to an effective exchange field, simple rate 
equation arguments are not sufficient. 

The main objective of the present paper is thus a sys- 
tematic analysis of charge and spin transport through a 
SMM. This has been achieved by employing the real-time 
diagrammatic technique^ which enables accurate study 
of transport properties in the full weak coupling regime. 
In particular, including the first- and second-order self- 
energy diagrams, we calculate the current, tunnel magne- 
toresistance (TMR) and shot noise in the presence of se- 
quential tunneling, cotunneling and cotunneling-assisted 
sequential tunneling processes. We show that the second- 
order processes determine transport in the Coulomb 
blockade regime, leading for instance to enhanced tun- 
nel magnetoresistance effect as compared to the value 
based on the Julliere model^ and to super-Poissonian 
shot noise due to bunching of inelastic cotunneling pro- 
cesses through the molecule. In addition, we also discuss 
the effects due to external magnetic field as well as the 
role of strength and type of exchange interaction between 
the molecule's spin and conduction electrons. 

The paper is organized as follows. In section II we de- 
scribe the model of a single-molecule magnet coupled to 
ferromagnetic leads. The real-time diagrammatic tech- 
nique used in calculations is presented briefly in section 
III. Section IV is devoted to numerical results and their 
discussion. In particular, the conductance, tunnel mag- 
netoresistance and shot noise in the linear and nonlinear 
response regimes are analyzed in subsections A and B, 
respectively. The dependence of transport properties on 
the strength of exchange coupling is discussed in subsec- 
tion C, while the effects of longitudinal external magnetic 
field are considered in subsection D. Furthermore, we also 
briefly discuss transport characteristics in the case when 
the exchange coupling is antiferromagnetic, subsection E. 
Finally, conclusions are given in section V. 



II. DESCRIPTION OF MODEL 

In this paper we consider a model SMM which is 
attached to two metallic ferromagnetic electrodes, see 
Fig. [T] The molecule is assumed to be weakly coupled 
to the leads, whose magnetizations form a collinear con- 
figuration, either parallel or antiparallel. The limit of 
strong coupling, where interesting phenomena such as the 
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Figure 1: (color online) Schematic representation of the sys- 
tem under consideration. The system consists of a SMM 
weakly coupled to two ferromagnetic electrodes with the 
collinear configuration of their magnetic moments, i.e. either 
parallel or antiparallel. Due to symmetrically applied bias 
voltage V — (fiL — fJ,R)/e, where ^l(r) denotes the electro- 
chemical potential of the left (right) lead, the LUMO level is 
independent of V. Position of the LUMO level, however, can 
be tuned by the gate voltage V g . 



Kondo effect22r— can be observed, is not considered here. 

Electronic transport through the molecule is assumed 
to take place only via the lowest unoccupied molecular or- 
bital (LUMO) of the SMM, which is coupled to the inter- 
nal magnetic core of the molecule via exchange interac- 
tion. Moreover, we also neglect all other unoccupied lev- 
els which are assumed to be well above the LUMO level 
and therefore cannot take part in transport for voltages 
of interest^ 8 - Furthermore, following previous theoretical 
studies^r— we restrict our considerations to the case of 
molecules with vanishingly small transverse anisotropy. 

Taking the above into account, a SMM coupled to ex- 
ternal leads can be described by Hamiltonian of the gen- 
eral form 



n = n 



SMM 



H 



leads 



(i) 



The first term on the right hand side describes the SMM 
and is assumed in the following form: 



HsMM = ~ 



s ^ 



[ D + Y, Dic - c - +D 2 c \ c A c i 

a 

g(M B (S z + s z )H z . (2) 



The first line of Eq. @ accounts for the uniaxial mag- 
netic anisotropy of a SMM, characterized by the uniax- 
ial anisotropy constant D of a free-standing (neutral) 
molecule. When a bias voltage is applied, the LUMO 
level can be charged with up to two electrons, which in 
turn can affect the magnitude of the uniaxial anisotropy. 
The relevant corrections are included by the constants 
D\ and D^. Moreover, S z denotes the z component of 
the internal (core) spin operator S, whereas c\ (c CT ) is 
the creation (annihilation) operator of an electron in the 
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LUMO level. We note that the Hamiltonian (2) is ap- 
plicable to situations where electronic structure of the 
molecules's magnetic core is not changed by adding one 
or two electrons to the LUMO level, except for modifica- 
tion of the anisotropy constants. 

The second line of the Hamiltonian Hsmm describes 
the LUMO level of energy e, with U being the Coulomb 
energy of two electrons of opposite spins that can oc- 
cupy this level. Although the position of the LUMO level 
can be modified by the gate voltage V g , it remains inde- 
pendent of the symmetrically applied bias voltage V. An 
important term for the present discussion is the last one, 
given explicitly by 



Js ■ S 



J 



J 



c[c t S + 



J 



4 c t 



Sz, (3) 



which stands for exchange coupling between the magnetic 
core of a SMM, represented by the spin S, and electrons 
in the LUMO level, described by the local spin opera- 
tor s = | J2aa' cJ-CTcto-'Cq-. where <r is the vector of Pauli 
matrices. This interaction can be either of ferromagnetic 
(J > 0) or antiferromagnetic (J < 0) type. Finally, the 
last term of Hsmm describes the Zeeman splitting associ- 
ated with the magnetic field applied along the easy axis 
of the molecule, where g stands for the Lande factor, and 
fiB is the Bohr magneton. 

In general, the molecular Hamiltonian "Hsmm is not 
diagonal, except for the case of a free-standing (un- 
charged) uniaxial SMM. It has been show n 25 i 26 i 28 that 
for the molecules with no transverse anisotropy, "Hsmm 
commutes with the zth component (S%) of the total spin 
S t = S + s, hence allowing us to analytically diagonal- 
ize it in the basis represented by the eigenvalues m of 
S* and the corresponding occupation number n of the 
LUMO level. In a general case, on the other hand, the 
problem can be dealt with numerically by performing a 
unitary transformation U^HsmmU = Hsmm to a new ba- 
sis in which Hsmm is diagonal. Consequently, we obtain 
the set of relevant eigenvectors \x) and the corresponding 
eigenvalues e x satisfying HsmmIx) = £ x\x)- 

The second term of Eq. describes ferromagnetic 
electrodes, and the gth electrode (q = L, R) is character- 
ized by noninteracting itinerant electrons with the dis- 
persion relation e^ CT , where k denotes a wave vector and 
a is the electron's spin. As a result, the lead Hamiltonian 
can be written as 



Hi, 



cads 



1 k,<r 



(4) 



where d£ (o^o-) ^ s ^ ne creation (annihilation) operator 
for an electron in the qth electrode. The degree of spin 
polarization of the ferromagnetic lead q can be described 
by the parameter P q , P q = (D\ - D c L)/{D q + +D q _), with 
D q ± denoting the density of states for majority (upper 
sign) and minority (lower sign) electrons at the Fermi 
level in the lead q. 



Finally, the last term Ht of the total Hamiltonian ((T|) 
describes tunneling processes between the molecule and 
the leads, and it is given by 



Hn 



E£[ r « 



T*A n q 



(5) 



with T q denoting the tunnel matrix element between the 
molecule and the qth lead. Due to the tunneling pro- 
cesses, the LUMO level of the molecule acquires a fi- 
nite spin-dependent width, T a = 5^_r£, where = 
2n\Tq\ 2 D q r . The parameters T q ± can be also expressed 
in terms of the spin polarization P q of the lead q as 
T q ± = T q (l±P q ) for spin-majority (upper sign) and spin- 
minority (lower sign) electrons, where T q — (T q ^+T q _)/2. 
In the following these parameters will be used to describe 
the strength of coupling between the LUMO level and the 
leads. Unless stated otherwise, the couplings are assumed 
to be symmetric, Tl = Tr = T/2. 



III. METHOD OF CALCULATIONS 

Among different available methods, only a few en- 
able us to analyze spin-dependent transport of the con- 
sidered system in both the sequential and Coulomb 
blockade regimes within one fully consistent theoreti- 
cal approach^ In particular, here, we employ the real- 
time diagrammatic techniquej 37 i 41 ' 44 ~ — which has al- 
ready proven its reliability and versatility in studying 
transport properties of various nanoscopic systems. 

The basic idea of this technique relies on a system- 
atic perturbation expansion of the reduced density ma- 
trix of the system under discussion and the operators of 
interest with respect to the coupling strength T between 
the LUMO level and the leads. All quantities, such as 
the current /, differential conductance G and the (zero- 
frequency) current noise S are essentially determined by 
the nonequilibrium time evolution of the reduced density 
matrix for the molecule's degrees of freedom. In the case 
considered in this paper, the density matrix has only di- 
agonal matrix elements, p x (t), which correspond to prob- 
ability of finding the molecule in state \x) at time t. Fol- 
lowing the matrix notation introduced by Thielmann et 
al.j& the vector p(t) of the probabilities is given by the 
relatior^Lii 



P (t) = n(i,t )p(t ). 



(6) 



where TL(t, to) is the propagator matrix whose elements, 
Tl x ' x (t, to), describe the time evolution of the system that 
propagates from a state \x) at time to to a state \x') at 
time t, and p(to) is a vector representing the distribu- 
tion of initial probabilities. In principle, the whole dy- 
namics of the system is governed by the time evolution 
of the reduced density matrix. Furthermore, this time 
evolution can be schematically depicted as a sequence 
of irreducible diagrams on the Keldysh contour^ which 
after summing up correspond to irreducible self-energy 
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blocks W x , x (i',i)4a The self-energy matrix W (*',*) is 
therefore one of the central quantities of the real-time di- 
agrammatic technique, as its elements W x ' x (£',£) can be 
interpreted as generalized transition rates between two 
arbitrary molecular states: \x) at time t and \x') at time 
t' . Consequently, the Dyson equation for the propagator 
is obtained in the for m 37 ' 41 ' 44 ' 45 

U(t,t ) = l+ [ dt 2 f diiW(f a ,*i)n(ti,to)- (7) 

By multiplying Eq. ([7|) from the right hand side with 
p(to)> an d differentiating it with respect to time t, one 
gets the general kinetic equation for the probability vec- 
tor p(i), 



(8) 



In the limit of stationary transport the aforemen- 
tioned formula reduces to the steady state master-like 
equation^iiiiii^ 



K), 



= r<5 



XX 



(9) 



where p st = lim t _ i . 00 p(t) = lim t0 _ i ._ oo p(0) is the station- 
ary probability vector, independent of initial distribution. 
On the other hand, W denotes the Laplace transform of 
the self-energy matrix W(t',f), whose one arbitrary row 
Xo has been replaced with (T, . . . , T) to include the nor- 
malization condition for the probabilities YlvPx = !■ 
Knowing the probabilities, the electric current flowing 
through the system can be calculated from the formula^ 



1 = 7rr Tr 



wV 



(10) 



where the matrix W i denotes the self-energy matrix in 
which one internal vertex originating from the expansion 
of tunneling Hamiltonian Ht has been substituted with 
an external vertex for the current operator. 

In order to calculate the transport quantities in both 
the deep Coulomb blockade and the sequential tunnel- 
ing regime in each order in tunneling processes, we per- 
form the perturbation expansion in T adopting the so- 
called crossover perturbation scheme^ i.e. we expand 
the self-energy matrices, W = Y^Li W n ) and W 7 = 
X]^Li W 1 ^. Here, the first order of expansion (n = 1) 
corresponds to sequential tunneling processes, while the 
second-order contribution (n = 2) is associated with co- 
tunneling processes. In the present calculations we take 
into account both the first- and second-order diagrams, 
which allows us to resolve the transport properties in the 
full weak coupling regime, i.e. in the cotunneling as well 
as in the sequential tunneling regimes. Furthermore, by 
considering the n = 1 and n = 2 terms of the expan- 
sion, we systematically include the effects of LUMO level 
renormalization, cotunneling-assisted sequential tunnel- 
ing, as well as effects associated with an exchange field 



exerted by ferromagnetic leads on the molecule^ 7 - — For 
n < 2, the stationary probabilities can be found from 
Eq. ©, with W = W< 4) +W< 2 ). On the other hand, the 
current is explicitly given by Eq. (|10j) where one has to 
take W J = W 7(1) + W /(2) . The key problem is now the 
somewhat lengthy but straightforward calculation of the 
respective self-energy matrices, which can be done using 
the corresponding diagrammatic rules . 37 ! 41 ' 44 ' 45 ' 47 An ex- 
ample of explicit formula for a second-order self-energy 
between arbitrary states |x) and \x'} can be found in 
Ref. 

With recent progress in detection of ultra-small signals, 
it has become clear that the information about the sys- 
tem transport properties can also be extracted from the 
measurement of current noise.— In fact, the shot noise 
contains information about various correlations, coupling 
strengths, effective charges, etc., which is sometimes un- 
accessible just from measurements of electric current. 
Therefore, to make the analysis more self-contained, in 
this paper we will also calculate and discuss the zero- 
frequency shot noise. The shot noise is usually defined 
as the correlation function of the current operators, and 
its Fourier transform in the limit of low frequencies is 
given by 5 - 1 S = 2/° <xj dt[(/(t)/(0) + /(0)/(t) - 2{/) 2 )]. 
For \eV\ > k^T, the current noise is dominated by fluctu- 
ations associated with the discrete nature of charge (shot 
noise), while for low bias voltages, the thermal noise dom- 
inates.— The general formula for the current noise within 
the language of real-time diagrammatic technique can be 
found in Ref. (H|. 



IV. NUMERICAL RESULTS AND DISCUSSION 

In this section we present and discuss numerical results 
on charge current, differential conductance, shot noise 
(Fano factor) and tunnel magnetoresistance (TMR) in 
the linear and nonlinear response regimes. The Fano fac- 
tor 



F 



2e J 



(11) 



describes deviation of the current noise from its Poisso- 
nian value, Sp = 2e\I\, which is characteristic of uncor- 
rected in time tunneling processes. On the other hand, 
the TMR is defined asi^ 7 -^ 



TMR = 



Ip - I At 
Iap 



(12) 



where Ip (Iap) is the current flowing through the system 
in the parallel (antiparallel) magnetic configuration at a 
constant bias voltage V. The TMR describes a change of 
transport properties when magnetic configuration of the 
device varies from antiparallel to parallel alignment - the 
conductance is usually larger in the parallel configuration 
and smaller in the antiparallel one, although opposite 
situation is also possible. 
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Figure 2: (color online) The total (first plus second order) 
differential conductance in the parallel (a) and antiparallel 
(b) configurations for the parameters: S = 2, J = 0.2 meV, 
D = 0.05 meV, £>i = -0.005 meV, D 2 = 0.002 meV, U = 1 
meV, fc B T = 0.04 meV, Pl = Pr = 0.5, and T = 0.002 meV. 



Numerical results have been obtained for a hypothet- 
ical SMM characterized by the spin number S = 2 and 
strong uniaxial magnetic anisotropy. However, we note 
that although in the following we assume S = 2, our con- 
siderations are still quite general and qualitatively valid 
for molecules with larger spin numbers. In fact, the choice 
of low molecule's spin allows us to perform a detailed 
analysis of various molecular states mediating the first 
and second-order tunneling processes. A large number of 
molecular states for S » 1 would make the discussion 
rather obscure. Apart from this, we assume a symmet- 
rical coupling of the molecule to the two external leads 
(Pl = Pr = P) and ferromagnetic exchange coupling be- 
tween the molecule's magnetic core and electrons in the 
LUMO level. Later on, however, we will relax the lat- 
ter restriction and consider the situation where the ex- 
change coupling is antiferromagnetic. For clarity reasons, 
we disregard the effects due to the negative sign of elec- 
tron charge, i.e. assume that charge current and particle 
(electron) current flow in the same direction (e > 0). 

We start from some basic transport characteristics of 
the system under consideration. In Fig. [5] we show the 
differential conductance in the parallel and antiparallel 
configurations as a function of the bias voltage and posi- 
tion of the LUMO level. The latter can be experimentally 
changed by sweeping the gate voltage. The density plot 
of the conductance displays the well-known Coulomb di- 
amond pattern. The average charge accumulated in the 
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Figure 3: (color online) Density plot of the total (first plus sec- 
ond order) TMR (a) and TMR calculated in the sequential 
tunneling approximation (b) plotted in the same scale and 
for the same parameters as in Fig. [21 The sequential TMR 
is smaller than the total TMR. The dashed lines are only a 
guide for eyes, and they represent positions of the main con- 
ductance peaks, Fig. [21 separating thus regions corresponding 
to different occupation states of the LUMO level. 



LUMO level is Q = ~}2, x n(x)p x (in the units of e), where 
n(x) = 0, 1, 2 denotes the number of additional electrons 
on the molecule in the state |x). When lowering energy of 
the LUMO level, the latter becomes consecutively occu- 
pied with electrons. This leads to two peaks in the linear 
conductance, separated approximately by U , which cor- 
respond to single and double occupancy, respectively, see 
Fig. H for V = 0. 

Furthermore, in the nonlinear response regime, the dif- 
ferential conductance shows additional lines due to tun- 
neling through excited states of the molecule. These fea- 
tures are visible in both magnetic configurations. On the 
other hand, the hallmark of spin-depended tunneling is 
the difference in magnitude of the conductance in paral- 
lel and antiparallel configurations - the conductance in 
the parallel configuration is generally larger than in the 
antiparallel one, see Fig. [5J This difference is due to spin 
asymmetry of tunneling processes, which leads to sup- 
pression of the conductance when configuration changes 
from parallel to antiparallel one. 

The density plot of the TMR corresponding to Fig. [5] 
is shown in Fig. [3]^ a). As one can note, the magnitude 
of TMR strongly depends on the transport regime. More 
precisely, TMR can range from approximately TMR w 
P 2 /(l - P 2 ) = 1/3 (for P = 0.5), which is characteris- 
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Figure 4: (color online) The average value of the zih compo- 
nent of the molecules's total spin (Si) for the parallel (a) and 
antiparallel (b) magnetic configurations. All parameters as in 
Fig.0 



tic of sequential tunneling regime where all states of the 
LUMO level are active in transport^ to roughly twice 
the value resulting from the Julliere model^ TMR 
TMR JulL = 4P 2 /(l - P 2 ) = 4/3, which can be observed 
in the nonlinear response regime of the Coulomb block- 
ade diamond (Q = 1), see Fig. [3Ja). For comparison, in 
Fig. GUb) we display the TMR calculated using only the 
sequential tunneling processes. One can see that the first- 
order TMR is generally smaller than the total (first plus 
second order) TMR. Furthermore, it is also clear that the 
second-order tunneling processes modify TMR mainly in 
the Coulomb blockade regime (Q = 1) as well as in the co- 
tunneling regimes where the LUMO level is either empty 
(Q = 0) or doubly (Q = 2) occupied. On the other hand, 
out of the cotunneling regime, the sequential processes 
dominate transport and the role of second-order tunnel- 
ing is relatively small. As a consequence, the two results 
become then comparable in these regions, see Fig. [3Ja) 
and Fig. |3][b). 

Spin-dependent transport through a SMM has a sig- 
nificant impact on its magnetic state. In Fig. U we show 
the average value of the molecule's spin zth component in 
the stationary state, (Sf), calculated as a function of the 
bias voltage V and energy of the LUMO level e. In the 
antiparallel magnetic configuration, Fig. Efb), the orien- 
tation of the molecule's spin is straightforwardly related 
to the bias voltage, and for V > the spin is aligned along 
the easy axis +z, whereas for V < it is aligned along 
the — z axis. Note, that in the regions corresponding to 
Q = and Q = 2 the spin is equal to that of magnetic 
core, while for Q = 1 it also includes the contribution 
from an electron in the LUMO level. By contrast, in the 
parallel configuration, Fig. @Ja), the value of (Sf) in the 
stationary state can be both positive and negative for 
each sign of the bias voltage, and it varies in a rather 
limited range close to zero. Moreover, (Sf) in the paral- 
lel (antiparallel) magnetic configuration is an even (odd) 
function of the bias voltage V. 

To account for the transport properties in different 
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Figure 5: (color online) Energy spectrum of the molecule un- 
der consideration (relevant parameters given in the caption 
of Fig. |5J for e = -0.5 meV (a,c) and e = 0.75 meV (b,d) in 
the case of ferromagnetic (a)-(b) and antiferromagnetic (c)- 
(d) coupling between the SMM's core spin and the spin of 
electrons in the LUMO level. The dashed line represents the 
Fermi level of the leads when no external voltage bias is ap- 
plied (V — 0). Different sets of molecular states correspond 
to different values of the SMM's total spin St, and/or the oc- 
cupation number of the LUMO level: |2; 0, m) (•), 1 5/2; I, m) 
(♦), |3/2;l,m) (A), and |2;2,m) (■). Note that in (a) and 
(c) the degeneracy between states |2; 0, m) and |2; 2, m) takes 
place only for m — 0. 



regimes, especially of TMR and shot noise, in the fol- 
lowing we present and discuss the gate and bias voltage 
dependence corresponding to various cross-sections of the 
relevant density plots mentioned above. More specifically, 
we will first consider transport properties in the linear 
response regime (Sec. IIV A|) . and then transport in the 
nonlinear regime (Sec. IIVB|) . In addition, whenever ad- 
visable and possible, we will also compare and relate our 
findings to existing results on quantum dot systems. At 
this point, it is however worth emphasizing that the prob- 
lem of electron transport through a SMM is much more 
complex and physically richer than in the case of single 
quantum dots^ This is because now the transfer of elec- 
trons occurs through many different many-body states of 
the coupled LUMO level and molecule's magnetic core, 
see Eq. ©. 

Since transport properties of a system are determined 
by its energy spectrum, it is instructive to analyze it in 
more detail. For molecules with only uniaxial anisotropy 
considered in this paper, the molecule's Hamiltonian 
Hsmm can be diagonalized analytically (the relevant for- 
mulas can be found in Ref. [28[ , from where the notation 
for molecular states has also been adopted) . Energy spec- 
trum of the molecule under consideration is presented 
in Fig. \5\ for two different values of the LUMO level 
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energy e and two values of the coupling parameter J. 
Each molecular state \St; n, m) is labelled by the total 
spin number St, the occupation number n of the LUMO 
level, and the eigenvalue m of the zth component of the 
molecule's total spin, S^ = S z + ^(cjc-f- — cjcjj, where 
the second term stands for the contribution from elec- 
trons in the LUMO level. The change of the LUMO level 
energy leads to the change in the energetic position of the 
spin-multiplets |5/2;l,m), |3/2;l,m) and |2;2,m) with 
respect to 2;0, m). The latter multiplet corresponds to 
uncharged molecule and therefore is independent of e, see 
Fig. El 



A. Transport in the linear response regime 

In this subsection we will focus on transport in the 
linear response regime. As we have already mentioned 
above, conductance in the linear response regime (see 
Fig. 2 for V = 0), displays two resonance peaks separated 
approximately by U. For J > and D(2S - 1) > k B T, 
one can assume that the molecule is in the spin states 
of lowest energy. The position of the conductance peaks 
(resonances) corresponds then to e = eoi, 



£01 



JS 
2 



'■^+D 1 S 2 



9Vb\H z \ 



(13) 



for the transition from zero to single occupancy of the 
LUMO level, and to e = £12, 



£12 



- J -l-u 

2 



(Da + D 2 )S 2 - 



9Vb\H z \ 



(14) 



for the transition from single to double occupancy. It is 
worth noting that the above expressions may be useful 
for estimating the coupling constant J from transport 
measurements. Moreover, from the above formulas one 
can conclude that the middle of the Coulomb blockade 
(Q=l in Fig. [J) regime corresponds to e = £ m , with 



U 
' 2" 



2Di 



D2 c2 



(15) 



which for the parameters assumed in calculations gives 
e m = —0.516 meV. Interestingly, e m is independent of 
the exchange coupling J, anisotropy constant D, and ex- 
ternal magnetic field H z , but it depends on the Coulomb 
interaction U, corrections D\ and D2 to the anisotropy 
due to finite occupation of the LUMO level, and the 
molecule's spin number S. In fact, owing to finite con- 
stants D\ and D2, the particle- hole symmetry is broken, 
which manifests itself in an asymmetric behavior of trans- 
port properties, as will be shown below. 

Figure IBfa) shows the total TMR in the linear re- 
sponse regime, where for comparison TMR in the sequen- 
tial transport regime is also displayed (dash-dotted line) . 
Clearly, the results obtained within the sequential tun- 
neling approximation, which yield a constant TMR equal 
to P 2 /(l — P 2 ), are not sufficient as the total (first plus 
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Figure 6: (color online) (a) TMR in the linear response regime 
for the parameters as in Fig. [2] (solid line). The dot-dashed 
line shows the TMR calculated in the first-order approxima- 
tion, (b) Average value of the zth component of the molecule's 
total spin in the (P) parallel (solid line) and (AP) antiparallel 
(dashed line) magnetic configurations. The dotted lines in (a) 
and (b) correspond to the case of D\ — D2 = 0. 



second order) linear TMR displays a nontrivial depen- 
dence on the gate voltage. This behavior stems from the 
dependence of the second-order processes on the occupa- 
tion number of the LUMO level. 

Generally, cotunneling processes can be divided into 
two groups with respect to whether or not the molecule 
remains in its initial state after a cotunneling process, 
Fig. [3a)-(b). Although the cotunneling events do not 
change the charge state of the molecule, they can, how- 
ever, modify its spin state (inelastic cotunneling). More- 
over, the inelastic cotunneling processes can lead to mag- 
netic switching of the molecule's spin between two low- 
est energy states, as shown schematically in Fig. [7jb). 
We note that in addition to double-barrier cotunneling 
processes which transfer charge between two different 
electrodes, there are also single-barrier cotunneling pro- 
cesses, where an electron involved in the cotunneling pro- 
cess returns back to the same electrode. Although the 
latter processes do not contribute directly to the current 
flowing through the system, they can affect all the trans- 
port properties in an indirect way, by altering spin state 
of the molecule. 



1. Cotunneling regime with empty and doubly occupied 
LUMO level 

When the LUMO level is either empty (Q = 0) or fully 
occupied (Q = 2), the total TMR in the correspond- 
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Figure 7: (color online) (a)-(b) Schematic representation of 
the elastic and inelastic electron cotunneling processes. The 
two bottom panels show examples of inelastic cotunneling pro- 
cesses leading to increase of the zth component of the SMM's 
spin in the situation when the LUMO level is: (c) empty 
(Q = 0) and (d) doubly occupied (Q = 2). Note that V ->■ 
in the linear response regime, so we assume = I-Ir + + . 
When in the ground state the molecule is occupied by a sin- 
gle electron, Q = 1, inelastic processes can generally occur via 
two virtual states associated with empty and doubly occupied 
LUMO level. 



ing cotunneling regions is slightly larger than the Jul- 
liere valued TMR JulL = 2P 2 /(l - P 2 ) (TMR Ju11 ' = 2/3 
for P = 0.5), see Fig. EJa). Electron transport in these 
regions is primarily due to elastic cotunneling processes 
which change neither the electron spin in the LUMO level 
nor the spin of molecule's core, and thus are fully coher- 
ent. An example of such process is sketched in Fig.JTJa). 
The enhancement of TMR above the Julliere value is then 
associated with the exchange coupling of the LUMO level 
to the molecule's core spin, which additionally admits in- 
elastic cotunneling processes in these regions. In addition, 
the enhanced TMR may also result from the fact that by 
using the crossover perturbation scheme^ 7 - we also in- 
clude some effects associated with third-order processes, 
which may further increase the TMR. Moreover, unlike 
the case of a single quantum dot , 47 ' 53 the maximal values 
of TMR reached for Q = and Q = 2 do not necessarily 
have to be equal, see Fig. Eta). Below, we discuss these 
new features in more details. 

From the energy spectrum displayed in Fig. [5jb) fol- 
lows that the dominant elastic transfer of electrons be- 
tween the leads for Q = takes place via the follow- 
ing virtual transitions: |2;0, — 2) O 5/2; L, — 5/2) and 
|2;0, 2) 5/2; 1,5/2) [indicated with black arrows in 
Fig-E^b)]. In the parallel configuration, the former tran- 
sitions establish the transport channel for minority elec- 
trons, whereas the latter ones for majority electrons. The 



asymmetry between the occupation probabilities of the 
states |2; 0, -2) and |2; 0, 2) (with |2; 0, 2) being favored), 
which occurs due to inelastic cotunneling processes, gives 
rise to increased transport of majority electrons. On the 
other hand, there is no such asymmetry in the antiparal- 
lel configuration. This, in turn, leads to an enhancement 
of the TMR above the Julliere value, Fig. EJa). 



Similar arguments also hold for the case of Q = 2, 
where the molecular states |2; 2, m) correspond to double 
occupancy of the LUMO level. The main difference as 
compared to the situation discussed above is that now in 
a cotunneling process the electron first has to tunnel out 
of the LUMO level and then another electron can tunnel 
onto the molecule [Fig.[7td)]. Analysis similar to that for 
Q = shows that in the parallel configuration the in- 
elastic cotunneling processes result in lowering of the zth 
component of the SMM's spin, see Fig. [6jb) . Moreover, 
the asymmetry between the occupation probabilities of 
the states |2; 2, —2) and |2;2,2), where now |2; 2, —2) is 
favored, leads to increased elastic cotunneling of spin ma- 
jority electrons and therefore gives rise to enhanced TMR 
for Q = 2. 



Another interesting feature of TMR in the linear re- 
sponse regime, shown in Fig. Eta), is the difference in its 
magnitude in the cotunneling regions corresponding to 
Q = and Q = 2. This is contrary to the case of Ander- 
son model, where the linear TMR was found to be sym- 
metric with respect to the particle- hole symmetry point, 
e = — C//24 7 - In the case considered here, the situation 
is different due to coupling of the LUMO level to the 
molecule's spin, and also due to occupation dependent 
corrections to the anisotropy constant, see Eq. ([2]). These 
corrections reduce the uniaxial anisotropy of the molecule 
with increasing number of electrons in the LUMO level. 
As a result, the height of the energy barrier between 
the two lowest molecular spin states is also diminished 
for Q = 1 and Q = 2, and so are the energy gaps be- 
tween neighboring molecular states within the relevant 
spin multiplets. For this reason, the probability distri- 
bution of the molecular states for Q = 2 (and also for 
Q = 1) is more uniform than for Q = 0, see the solid line 
in Fig. Etb). Consequently, the value of TMR for Q = 2 
is smaller than for Q = 0. Thus, the observed asymmetry 
with respect to e = e m is due to the lack of particle-hole 
symmetry in the system when D\ and D2 are nonzero. 
However, if the influence of the LUMO level's occupa- 
tion on the anisotropy were negligible, D± pa D2 ~ 
(the states 2;0,m) and |2;2,m) in Fig. EKa) were then 
degenerate for every m), the symmetry with respect to 
£ = £m = —U/2 would be restored. This situation is pre- 
sented by the dotted curves in Fig. [6j which clearly show 
that the asymmetric behavior of TMR and (S^) is related 
to the corrections to anisotropy constants and the lack of 
particle-hole symmetry. 
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2. Cotunneling regime with singly occupied LUMO level 

Interestingly, in the Coulomb blockade regime, with 
one electron in the LUMO level (Q = 1), the TMR 
reaches local maxima close to the center of the Coulomb 
gap, and a shallow local minimum just in the middle, 
i.e. for e = s m . This behavior is different from that ob- 
served in single-level quantum dots, where linear TMR 
in the Coulomb blockade regime becomes suppressed and 
reaches a global minimum when e = — As in the 

case of Q = and Q = 2 discussed above, the origin of in- 
creased TMR for Q = 1 can be generally assigned to the 
modification of the probability distribution of molecular 
states due to inelastic cotunneling processes, Fig. EJb). 
In turn, the appearance of the local minimum in the cen- 
ter of the Q = 1 region is related to the fact that when 
e = e m , the virtual states for leading inelastic cotunnel- 
ing processes, which belong to spin multiplets |2;0,m) 
and 1 2; 2, m), become pairwise degenerate (in the present 
situation, |2;0, ±2) with |2;2,±2)). This means that in 
the parallel configuration cotunneling processes involv- 
ing empty and doubly occupied virtual states occur at 
equal rates. As a consequence, the average spin on the 
molecule tends to zero, see Fig. EJb), and TMR displays 
a local minimum for e = e m . 



3. Resonant tunneling regime 

For resonant energies, Eqs. (fT3 |) - (fT4|) . where the oc- 
cupancy Q of the molecule changes, the sequential tun- 
neling processes play a dominant role. This results in 
the reduction of TMR to approximately half of the Jul- 
liere valued see the boundaries between the hatched and 
non-hatched areas in Fig. [6] The rate of first-order tun- 
neling processes increases whenever the two neighboring 
charge states of the molecule become degenerate, pro- 
vided that the conditions |An| = 1 and |ASf | = 1/2 are 
simultaneously satisfied, where |An| and |AS*| describe 
change in the occupation and spin of the molecule. This 
means that for e = £oi ~ 0-18 meV the degeneration be- 
tween the empty and singly occupied states, |2; 0; ±2) and 
1 5/2 ; 1, ±5/2), is observed, whereas for e = £12 ~ —1.21 
meV the states with a single and two electrons on the 
LUMO level, |5/2;l,±5/2) and |2;2;±2), are degener- 
ate. Moreover, we also note that for T ss /c^T, TMR can 
be reduced further due to increased role of second-order 
processes giving rise to the renormalization of the LUMO 
levels 



B. Transport in the nonlinear response regime 

The influence of sequential tunneling on transport 
characteristics, as well as on magnetic state of the SMM, 
grows with increasing bias voltage. For voltages above 
the threshold for sequential tunneling, first-order pro- 
cesses determine transport and the influence of cotun- 
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Figure 8: (color online) Bias dependence of the current (a,e), 
differential conductance (b,f), Fano factor (d,h) in the parallel 
(solid lines) and antiparallel (dashed lines) configurations and 
TMR (c,g) for e = -0.5 meV (a)-(d) and e = 0.75 meV (e)- 
(h). The parameters are the same as in Fig. [2l and Iq = 
eT/h m 0.5 nA. The dotted lines show the results obtained 
taking into account only first-order tunneling processes. The 
effect of cotunneling is most pronounced in the TMR and 
Fano factor. 



neling is rather small. However, when applied voltage is 
below the threshold, sequential tunneling becomes ex- 
ponentially suppressed and second-order processes give 
the dominant contribution to the current, and need to 
be taken into account to get a proper physical picture. 
Figure [5] shows the bias dependence of the current, differ- 
ential conductance, TMR and Fano factor, calculated for 
e = —0.5 meV and e = 0.75 meV. The former case cor- 
responds to the situation where the LUMO level in equi- 
librium is singly occupied, Fig. [Ha) , while in the latter 
case it is empty, Fig.[5jb). One can see that cotunneling 
significantly modifies the first-order results in the block- 
ade regimes and this modification is most pronounced for 
TMR and shot noise. 



1. Transport characteristics in the case o/eoi > e > £12 

Consider first the case when in equilibrium the LUMO 
level is singly occupied (left panel of Fig. [5]) . At low tem- 
peratures and low voltages, the molecule with almost 
equal probabilities is in one of the two ground states 
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1 5/2 ; 1, ±5/2), Fig. H^a). When a small bias voltage is 
applied, some current flows due to cotunneling processes 
through virtual states of the system. If the bias volt- 
age exceeds threshold for sequential tunneling, the cur- 
rent significantly increases and becomes dominated by 
first-order processes, when electrons tunnel one-by-one 
through the molecule. 

Since elastic cotunneling in the antiparallel configu- 
ration occurs essentially through the minority-majority 
and majority- minority channels, whereas for parallel 
alignment through the majority-majority and minority- 
minority ones, one observes growth of TMR with in- 
creasing bias voltage, which reaches a local maximum 
just before the threshold for sequential tunneling. This 
is associated with nonequilibrium spin accumulation in 
the LUMO level for the antiparallel configuration, which 
leads to suppression of charge transport and thus to en- 
hanced TMR. Further increase of transport voltage re- 
sults in a decrease of TMR to approximately 1/3 (for 
Pl = Pr = 0.5), which is typical of the sequential tun- 
neling regime, when all molecular states actively partic- 
ipate in transport ^ 7 ' 53 In the parallel magnetic configu- 
ration all states are then equally populated, so that aver- 
age magnetic moment of the molecule vanishes, (S^)=0. 
This differs from the antiparallel case, in which only the 
states with large positive zth component of the SMM's 
spin have remarkable probabilities. Finally, we note that 
the slight shift between the peaks in differential conduc- 
tance corresponding to different magnetic configurations, 
see Fig.[8jb), is a consequence of nonequilibrium spin ac- 
cumulation in the LUMO level in the antiparallel con- 
figuration. Similar behavior has been observed in the 
case of transport through ferromagnetic single-electron 
transistors 44 

The Fano factor in the parallel (Fp) and antiparal- 
lel (Fap ) configurations is shown in Fig. [Sf d) . For low 
bias voltages, the shot noise is determined by thermal 
Johnson-Nyquist noise, which results in a divergency of 
the Fano factor for V —> (current tends to zero). 
When a finite bias voltage is applied to the system, the 
Fano factor in both magnetic configurations drops to the 
value close to unity, which indicates that transport oc- 
curs mainly due to elastic cotunneling processes. Such 
processes are stochastic and uncorrelated in time, so the 
shot noise is Poissonian. When bias voltage increases fur- 
ther, the shot noise is enhanced due to bunching of in- 
elastic cotunneling processes and reaches maximum just 
before threshold for sequential tunneling. At the thresh- 
old voltage, sequential tunneling processes begin to dom- 
inate transport and the noise becomes sub-Poissonian. 
This indicates that tunneling processes in the sequen- 
tial tunneling regime are correlated due to Coulomb cor- 
relation and Pauli principle, which generally gives rise 
to suppressed shot noise as compared to the Poissonian 
value. Furthermore, another feature clearly visible in the 
Coulomb blockade regime is the difference in Fano fac- 
tors for parallel and antiparallel magnetic configurations. 
More specifically, shot noise in the parallel configura- 



tion is larger than in the antiparallel one. This behavior 
is associated with the fact that transport in the paral- 
lel configuration occurs mainly through two competing 
majority-majority and minority-minority spin channels, 
which in turn increases fluctuations, thus Fp > Fap- 



2. Transport characteristics in the case of e > eoi 

Let us consider now the situation shown in the right 
panel of Fig. i.e. when the LUMO level of the molecule 
is empty at equilibrium, Fig.^b). The initial large value 
of TMR, whose origin was discussed above, drops sharply 
as the bias voltage approaches the threshold value for se- 
quential transport. In turn, the first pronounced peak 
in differential conductance appears when the following 
transitions become allowed: |2;0,±2) O |5/2;l,±5/2) 
[denoted by arrows is Fig. [Hb)]. It is important to note 
that, when a spin-multiplet enters the transport energy 
window, the first states that take part in transport are 
those with the largest \{S%)\ (lowest energy). Conse- 
quently, in the parallel magnetic configuration the sys- 
tem can be temporarily trapped in some molecular spin 
states of lower energy. For larger bias voltage, additional 
small peaks appear in the conductance for parallel con- 
figuration, and some of them are also visible in the an- 
tiparallel configuration. In general, these peaks are re- 
lated to transitions involving states from the multiplet 
|3/2;l,m): |2;0,±1) O 3/2; 1, ±3/2) (A), |2;0,±2) o 
|3/2;l,±3/2) (B) and 3/2; 1, ±3/2) «-> |2;2,±2) (C), 
respectively, see Fig. Of). In the parallel configuration 
all the three peaks are visible, whereas for antiparal- 
lel alignment only the peak B can be clearly resolved. 
Since in the antiparallel configuration tunneling processes 
tend to increase the zth. component of the SMM's total 
spin, the probability of finding the molecule in any of the 
spin states |2;0, m) differs significantly from zero only 
for m = 2. As a consequence, in the antiparallel config- 
uration most favorable transitions are those having the 
initial state |2;0, 2), and thus the peaks A and C are 
suppressed, see Fig.JHJf). 

Furthermore, as soon as all states within a certain spin- 
multiplet become energetically accessible, the probability 
of finding the molecule in each of these states becomes 
roughly equal. On the other hand, in the antiparallel con- 
figuration the system tends towards maximum value (for 
V > 0) of the zth component of SMM's spin. For these 
reasons, some regions of the increased TMR are present 
in Fig. Big). 

The corresponding Fano factor is shown in Fig. [Stff). 
At low bias, the Fano factor drops with increasing volt- 
age. However, its bias dependence is distinctively differ- 
ent in both magnetic configurations. In the antiparallel 
configuration, the Fano factor tends to unity, indicating 
that transport is due to uncorrelated tunneling events. In 
the parallel configuration, on the other hand, we observe 
large super-Poissonian shot noise. The increased current 
fluctuations result mainly from the interplay between 
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Figure 9: (color online) Density plot of TMR for J = 0.05 
meV (a) and J = 0.4 meV (b). The bottom panels display 
cross-sections of (a) and (b) for the constant energy of the 
molecule's LUMO level e = —0.5 meV (c) and the constant 
bias voltage V = mV (d). Other parameters are the same 
as in Fig. [2] 



different cotunneling processes and bunching of inelas- 
tic cotunneling. In addition, as mentioned previously, in 
the parallel configuration the molecule can be temporar- 
ily trapped in some molecular spin states of lower en- 
ergy, which also gives rise to super-Poissonian shot noise. 
When the bias voltage is increased above the threshold 
for sequential tunneling, the Fano factor becomes sup- 
pressed and the shot noise is generally sub-Poissonian. 
Finally, we also note that super-Poissonian shot noise 
in the cotunneling regime has already been observed in 
quantum dots and carbon nanotubeS ) 53 i 55 ~ — where the 
increased noise was associated with bunching of inelastic 
spin-flip cotunneling events. 



level. As follows from Fig. Ufa), behavior of TMR for 
small values of J starts bearing some resemblance to 
that of a single-level quantum dotj 47 i 53 Furthermore, in 
the linear response regime, Fig.^d), the enhanced TMR 
around the electron-hole symmetry point is no longer vis- 
ible, and instead a global minimum develops there. In 
fact, in the limit of J = one observes a simple quantum- 
dot-like transport behavior. 47 ! 53 

In the opposite limit of large J shown in Fig. [^b), the 
maxima in the total linear TMR are shifted away from 
the zero bias point. This is a consequence of increased 
energy gaps between the ground states 5/2; 1, ±5/2) 
and the nearest lying states satisfying An = 1 and 
\AS?\ = 1/2, i.e. |2;0(2),±2). Another interesting fea- 
ture of TMR visible in the linear response regime is the 
presence of additional two local minima around e = e m . 
see Fig. [9jd). Some precursors of these minima can be 
actually seen also in Fig. HJa) as two steep steps on both 
sides of the plot's central part. Generally, they stem from 
an uneven probability distribution of the molecular spin 
states with positive and negative zth component of the 
SMM's spin in the parallel magnetic configuration, see 
Fig. IHlIb). This in turn means that elastic cotunneling 
processes occur mainly through the minority-minority 
spin channel, so that transport is effectively suppressed. 
In the present situation, the minima are more distinct due 
to larger energy separation between the spin-multiplets 
|5/2;l,m) and |3/2;l,m). 

In the nonlinear response regime, on the other hand, 
the TMR exhibits a minimum at zero bias and starts 
increasing with the bias voltage to reach a maximum 
around the threshold for sequential tunneling. This is as- 
sociated with nonequilibrium spin accumulation in the 
LUMO level, which is present in the antiparallel con- 
figuration. We note that although the magnitude and 
position of the TMR maxima in the nonlinear response 
regime of the Coulomb blockade depend significantly on 
the exchange constant J, the general qualitative behav- 
ior of TMR is rather independent of J, see Figs.^c) and 
He). 



D. Transport in the presence of a longitudinal 
external magnetic field 



C. Dependence on exchange coupling strength 

Tunnel magnetoresistance may become significantly 
changed by altering the strength of ferromagnetic ex- 
change coupling between the LUMO level and the SMM's 
core spin, as shown in Fig. [9] With decreasing </, the en- 
ergy separation between the relevant molecule states cor- 
responding to the single occupancy of the LUMO level is 
also diminished (slanted squares and triangles in Fig. [5] 
start then approaching each other). This, in turn, leads to 
a reduction in size of the central diamond-shaped region, 
representing transport in the Coulomb blockade regime 
through the molecule with one electron in the LUMO 



Let us consider now the main effects due to a finite 
magnetic field applied to the system. When the field is 
along the easy axis of the molecule, its effects occur via 
modification of the energy of molecular spin states. On 
the other hand, when the field possesses also a transver- 
sal component, it leads to symmetry-breaking effects and 
the zth component of the SMM's total spin is no more 
a good quantum number^ If the magnetic field is addi- 
tionally time-dependent, one can expect the phenomenon 
of quantum tunneling of magnetization to occur . 28 ' 59 ! 60 
Since the primary focus of the present paper is on trans- 
port through SMMs with uniaxial anisotropy, in the fol- 
lowing we consider only a longitudinal magnetic field. 
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Figure 10: (color online) (a) Density plot of TMR in the case 
when the external magnetic field H z = 0.216 T (gusHz = 
0.025 meV) is applied along the easy axis of the molecule, (b) 
TMR in the linear response regime (solid line). For compar- 
ison, TMR in the absence of external magnetic field (dotted 
line in (b)) is also shown. The other parameters are the same 
as in Fig. [2] 
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Figure 12: (online color) The current (a,e), differential con- 
ductance (b,f), Fano factor (d,h) in the parallel (solid lines) 
and antiparallel (dashed lines) configurations, and the TMR 
(c,g) for e = -0.5 meV (a)-(d) and e = 0.75 meV (e)-(h) 
as a function of the bias voltage. An external magnetic field 
H z — 0.216 T is applied along the z-axis, while the other 
parameters are the same as in Fig. [2] 



Figure 11: (color online) Average value of the zth compo- 
nent of the total molecule's spin in the parallel (a) and an- 
tiparallel (b) magnetic configurations, when an external field 
H z = 0.216 T is applied along the z-axis. The other parame- 
ters are the same as in Fig. [2] 



Figure [TW a) shows the density plot of TMR for a mag- 
netic field applied along the easy axis of a SMM. Despite 
rather modest value of the field (for comparison, in the 
experiment on the Mni2 molecule attached to nonmag- 
netic metallic electrodes by Jo et at, the field of 8 T 
was used, Ref. (ll|), a drastic change in transport prop- 
erties of the system is observed [contrast Fig. [TUta) with 
Fig. OJa)]. First, the field breaks the symmetry with re- 
spect to the bias reversal. Second, it admits the situation 
when conductance in the antiparallel magnetic configura- 
tion is larger than in the parallel one (black regions cor- 
responding to negative TMR). Furthermore, in the par- 
allel configuration the average spin (S^) in the Coulomb 
blockade region can take large negative values, while in 
the absence of magnetic field the SMM's spin prefers ori- 
entation in the plane normal to the easy axis. This implies 
that for parallel alignment of leads' magnetizations, the 



molecule's spin tends to align antiparallel to the z-axis, 
Fig. Illf a). However, when the sequential tunneling pro- 
cesses are allowed, this tendency is generally reduced. In 
the antiparallel configuration, on the other hand, the be- 
havior of the average molecule's spin is similar to that 
for H z = 0, see Figs.QjJb) andg^b). 

In the linear response regime, a large change of TMR 
is observed when e is comparable to e m , i.e. in the mid- 
dle of the Coulomb blockade regime, see Fig. [TOf b) . This 
stems from the fact that at this point the dominating 
spin-dependent channel for transport due to cotunneling 
processes in the parallel magnetic configuration switches 
from the minority-minority channel (for s > e m ) to 
majority-majority one (for e < e m ). In the antiparallel 
configuration, on the other hand, the dominant channel 
is rather associated with majority-minority spin bands, 
irrespective of the position of the LUMO level. As a con- 
sequence, for e > e m the current in the parallel configura- 
tion is smaller than that in the antiparallel one, leading 
to negative TMR, whereas for e < e m the situation is 
opposite and one finds a large positive TMR effect, see 
Fig. Mb). 

The transport characteristics in the nonlinear response 
regime, and in the presence of external magnetic field, are 
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shown in Fig. [T^l where the left (right) panel corresponds 
to the situation where in the ground state the molecule is 
singly occupied (empty). The asymmetry with respect to 
the bias reversal is clearly visible, especially in the tunnel 
magnetoresistance, see Fig. IT27 c) and (g). Interestingly, 
this asymmetric behavior is mainly observed in the co- 
tunneling regime, as can be also seen in Fig. UOf aL This 
results from the fact that the degeneracy of the molecule's 
ground state is removed for H z ^ and the SMM be- 
comes polarized. In turn, transport in the cotunneling 
regime depends mainly on the system's ground state, 
which is the initial state for the cotunneling processes. As 
a consequence, in the parallel configuration the current 
is always mediated by electrons belonging to the same 
spin bands of the leads, whereas in the antiparallel con- 
figuration, the dominant transport channel is associated 
either with majority or minority electrons, depending on 
the direction of the current flow. Thus, the current in the 
antiparallel configuration becomes in general asymmetric 
with respect to the bias reversal, which gives rise to the 
associated asymmetric behavior of TMR. 

For voltages larger than the splitting due to the Zee- 
man term (g/iBH z = 0.025 meV), the inelastic cotunnel- 
ing processes start taking part in transport. The com- 
petition between the elastic and inelastic cotunneling 
leads in turn to large super-Poissonian shot noise, which 
in the parallel configuration is enhanced due to addi- 
tional fluctuations associated with cotunneling through 
majority-majority and minority-minority spin channels, 
see Fig. [T27 d) and (h). On the other hand, when the 
voltage exceeds threshold for sequential tunneling, more 
states take part in transport and the asymmetry with 
respect to the bias reversal is suppressed. The same ten- 
dency is observed in the shot noise, which in the sequen- 
tial tunneling regime becomes generally sub-Poissonian. 



E. Antiferromagnetic coupling between the LUMO 
level and SMM's core spin 

The numerical results presented up to now concerned 
the case of ferromagnetic coupling ( J > 0) between the 
LUMO level and the SMM's core spin. However, since 
the type of such an interaction generally depends on the 
SMM's internal structure, the exchange coupling can be 
also of antiferromagnetic type (J < 0). In this subsec- 
tion we discuss how the main transport properties of the 
system change when the exchange coupling becomes an- 
tiferromagnetic. 

First of all, we note that in the case of anti- 
ferromagnetic coupling between the LUMO level and 
molecule's core spin the formulas estimating the posi- 
tion of conductance resonances need some modification. 
Equations (fT3")) - (fbI)) were derived assuming the degen- 
eracy between the states 2;0, ±2) (|5/2; 1, ±5/2)) and 
5/2; 1, ±5/2) (|2;2,±2)). For J < 0, however, the con- 
dition has to be modified by changing 1 5/2; 1, ±5/2) into 
1 3/2 ; 1, ±3/2), where the upper signs apply for H z < 0, 



TMR 



-0.4 -0.2 0.2 0.4 0.6 




-2-1012 -2-1012 
Bias voltage V (mV) Bias voltage V (mV) 

Figure 13: (color online) (a) The total tunnel magnetoresis- 
tance in the case of antiferromagnetic coupling between the 
SMM's core spin and the spin in the LUMO level, calculated 
for J = —0.2 meV and other parameters as in Fig. [2] (b) Rep- 
resentative cross-sections of the density plot in (a) for several 
values of the LUMO level energy e. 

and the lower ones for H z > 0. The relevant equations 
take now the following form: 

s 01 = ^ + D 1 S 2 - As (16) 

for the transition from empty to singly occupied states, 
and 

ei 2 = -^-U + (D 1 +D 2 )S 2 + As (17) 

for the transition between singly and doubly occupied 
states, where 

2 2 

- \JdM{DM + I J|)^^ + ^(2s7l)*, (18) 

with £>W = D + Di. 

The most apparent new feature of the total TMR for 
J < 0, as shown in Fig. U3f aL is its negative value 
in the Coulomb blockade regime (Q = 1). The nega- 
tive TMR occurs in transport regimes where the max- 
imum of TMR was observed for J > 0, i.e. close to the 
threshold for sequential tunneling, see Fig. Such 
behavior of TMR originates from the fact that now spin- 
multiplets |5/2;l,m) and |3/2;l,m) exchange their po- 
sitions, Fig.[^c)-(d), so that the multiplet corresponding 
to smaller total spin of the molecule for antiferromag- 
netic coupling corresponds to lower energy. Consequently, 
in the Coulomb blockade the current flowing in the an- 
tiparallel configuration is larger than that in the parallel 
configuration, which gives rise to the negative TMR ef- 
fect. 

The linear response TMR is shown in Fig.[T2|a). Unlike 
the case of ferromagnetic coupling, the values of TMR for 
Q = and Q = 2 are smaller as compared to those in the 
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Figure 14: (color online) Tunnel magnetoresistance (a) and 
the zth component of the molecule's total spin (b) calculated 
in the linear response regime for the antiferromagnetic cou- 
pling of the SMM's core spin with the spin of the LUMO level 
(J = —0.2 meV and other parameters as in Fig. [2j. Dotted 
lines show the results obtained for the case of ferromagnetic 
exchange coupling, see Fig. [6] - in (b) the dotted line corre- 
sponds to the parallel configuration. 



case of transport through single-level quantum dots i i 53 
On the other hand, for Q = 1 the TMR can take values 
exceeding those found in the case of ferromagnetic ex- 
change coupling. For s > Soi, the equilibrium probability 
distribution of different molecular spin states 2; 0, m) be- 
comes changed owing to inelastic cotunneling processes, 
similarly as described in Sec. IIV A[ The key difference 
with respect to J > is that now dominating elastic 
cotunneling transitions for Q = are those with ini- 
tial states |2;0,±2) and virtual states |3/2;l,±3/2) [in- 
dicated by black arrows in Fig.[5Jcl)]. 



V. SUMMARY AND CONCLUSIONS 

We have systematically analyzed the transport proper- 
ties of a single-molecule magnet coupled to ferromagnetic 
leads in both sequential and cotunneling regimes. The 
transport processes in such a system occur due to tunnel- 
ing through the LUMO level which is exchange coupled 
to the molecule's core spin. By employing the real-time 
diagrammatic technique, we have calculated the current, 
tunnel magnetoresistance and shot noise in both the lin- 
ear and nonlinear response regimes. The results show that 
the inclusion of second-order processes is crucial for a 
proper description of transport characteristics. 

Assuming the ferromagnetic coupling between the 
LUMO level and the molecule's core spin, we have shown 



that TMR in the Coulomb blockade regime can be en- 
hanced above the Julliere value. This enhancement is as- 
sociated with nonequilibrium spin accumulation in the 
molecule. Moreover, we have found an asymmetric behav- 
ior of the linear response TMR with respect to the middle 
of the Coulomb blockade regime, and its strong depen- 
dence on the number of electrons in the LUMO level. The 
asymmetry is associated with corrections to anisotropy 
constant due to a nonzero occupation of the molecule, 
which breaks the particle-hole symmetry in the system. 
In addition, we have shown that the competition between 
the elastic and inelastic second-order processes leads to 
large super-Poissonian shot noise. The shot noise is fur- 
ther enhanced in the parallel configuration due to addi- 
tional fluctuations associated with majority-majority and 
minority-minority spin channels for electronic transport. 
On the other hand, for bias voltages above the threshold 
for sequential tunneling, the shot noise becomes gener- 
ally sub-Poissonian, indicating the role of correlations in 
sequential transport. 

We have also discussed how transport properties de- 
pend on the strength of the exchange coupling J between 
the LUMO level and the molecule's core spin. When the 
exchange coupling is relatively weak, the transport be- 
havior of the system resembles that of single-level quan- 
tum dots, whereas with increasing exchange constant, 
the transport characteristics change in a nontrivial way 
and become distinctively different from those of quantum 
dots. In addition, it turned out that the position of max- 
ima of TMR in the Coulomb blockade depend linearly on 
the strength of the exchange coupling. This may be use- 
ful in determining the magnitude of exchange constant 
experimentally. 

Furthermore, we have studied the effects of external 
magnetic field and shown that current flowing through 
the SMM becomes then asymmetric with respect to 
the bias reversal. We have found a strong dependence 
of TMR on the number of electrons occupying the 
LUMO level. When the LUMO level is empty, the TMR 
may become negative, while for doubly occupied LUMO 
level tunnel magnetoresistance is much enhanced. Fi- 
nally, we have also discussed how transport properties 
change when the coupling between the LUMO level and 
molecule's core becomes antiferromagnetic. In that case 
we predict a large negative TMR effect in the Coulomb 
blockade regime, exactly where for ferromagnetic cou- 
pling an enhanced TMR was observed. Thus, the sign 
of TMR may provide an information on the type of ex- 
change interaction, which may be of assistance for future 
experiments. 

To conclude, we note that although the numerical re- 
sults presented in this paper concern SMMs coupled to 
ferromagnetic leads, most of the qualitative results are 
applicable also to SMMs coupled to nonmagnetic elec- 
trodes. Apart from this, we note that the model we have 
studied also corresponds to systems consisting of a single- 
level quantum dot exchange-coupled to a spin S. In fact, 
very recently a similar device built of a quantum dot 
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coupled through spin exchange interaction to metallic is- 
land have been implemented to experimentally access the 
quantum critical point between the Fermi liquid and non- 
Fermi liquid regimes^ 
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